Pseudogap induced by short-range spin correlations in a doped Mott insulator 
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We study the evolution of a Mott-Hubbard insulator into a correlated metal upon doping in the 
two-dimensional Hubbard model using the Cellular Dynamical Mean Field Theory. Short-range 
spin correlations create two additional bands apart from the familiar Hubbard bands in the spectral 
function. Even a tiny doping into this insulator causes a jump of the Fermi energy to one of these 
additional bands and an immediate momentum dependent suppression of the spectral weight at 
this Fermi energy. The pseudogap is closely tied to the existence of these bands. This suggests a 
strong-coupling mechanism that arises from short-range spin correlations and large scattering rates 
for the pseudogap phenomenon seen in several cuprates. 

PACS numbers: 71.10.Fd, 71.27.+a, 71.30.+h, 71.10.-w 



The issue of the origin of the pseudogap phenomenon 
observed in underdoped cuprates lies at the center of any 
theoretical explanation for high temperature supercon- 
ductivity in the cuprates and is one of the most challeng- 
ing questions in condensed matter physics. The suppres- 
sion of low energy spectral weight in the normal state of 
these materials has been observed through various exper- 
imental probes . In spite of many theoretical works to 
explain the observed anomalies, there is no consensus at 
present. The lack of controlled approximations to deal 
with the strong coupling physics and low dimensionality 
inherent to these systems continues to pose major stum- 
bling blocks towards a complete theoretical understand- 
ing. Since the parent compounds of the cuprates are 
Mott-Hubbard insulators, an understanding of such an 
insulator and its evolution into a correlated metal upon 
doping is crucial. 

In this paper we study the two-dimensional Hubbard 
model on a square lattice at and near half-filling with Cel- 
lular Dynamical Mean-Field Theory (CDMFT) 2]. The 
CDMFT method is a natural generalization of the single 
site DMFT Q to incorporate short-range spatial corre- 
lations. Since at and near half-filling short-range spin 
correlations are dominant at low energy, this method is 
expected to describe additional features caused by spin 
degrees of freedom in the single-particle spectrum. The 
CDMFT has already passed several tests against ex- 
act results obtained by the Bethe Ansatz and the Den- 
sity Matrix Renormalization Group (DMRG) techniques 
in one dimension, where the CDMFT scheme is expected 
to be in the worst case scenario. Long-range order in- 
volving several lattice sites such as ci-wave superconduc- 
tivity can be also described in CDMFTjj. Several other 
cluster schemes have been proposed 0, UH, El H3, includ- 
ing Dynamical Cluster Approximation (DCA) pfl, Clus- 
ter Perturbation Theory (CPT) Q and its variational 
extension (V-CPT) The variational principle used in 



the last scheme allows one to consider CPT, V-CPT, and 
CDMFT within a unified framework. 

In the CDMFT construction 0,0 the infinite lattice is 
tiled with identical clusters of size N c . In an effective ac- 
tion description, the degrees of freedom in a single cluster 
are treated exactly, while the remaining ones are replaced 
by a bath of non- interacting electrons, which hybridizes 
with the cluster degrees of freedom. For practical pur- 
poses, it is useful to view this cluster action as arising 
from a cluster-bath Hamiltonian of the form 

(liv),a A 1 
"F ^ y £m(T®rnaQ J 7nt7 "F ^ ^ Vm^icr (^rna^f 1 ^ H.C.) .(1) 

Here the indices /Lt, v — 1, • • • ,N C label sites within the 
cluster, to = 1, • • • , Nb with JV& representing the num- 
ber of bath degrees of freedom, and c^ a and a ma annihi- 
late electrons on the cluster and the bath, respectively. 
ta V is the hopping matrix within the cluster, e ma is the 
bath energy and V mtl(J is the bath-cluster hybridization 
matrix. Let us adopt a matrix notation (cluster indices 
suppressed) for the hopping matrix t, the cluster Green 
function G c {iuj), its nonintcracting counterpart Gq{iuj) 
and the cluster self-energy E = Gq 1 — G^ 1 . We start 
with an initial guess for the bath parameters e ma and 
V mtia which determines a starting Go(iui). With this 
guess the cluster Green function is calculated by solv- 
ing the cluster-bath Hamiltonian (Eq. ^| . To close the 
self-consistency loop we obtain a new G' (iuj) using 
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fi — t(k) — T,(iu> r , 
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The integral comes from projecting on the cluster the 
Green function of an infinite lattice formed of identical 
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clusters located at positions Ri . The self-energy of that 
Green function is equal to that in a single cluster. In real 
space it vanishes for sites that do not belong to the same 
cluster. However, hopping between clusters is allowed in 
the same way as on the original infinite lattice. Hence, 
the hopping matrix has a dependence on superlattice vec- 
tors Ri — Rj and correspondingly on wave vectors k that 
diagonalize that part of the problem. Now to invert the 
relation between this G' (iu>) and the bath parameters for 
the next iteration we minimize 
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The lattice Green function G(k, iuj n ) is obtained from 
G(k,iuj n ) = -jTT^t 



ik-(f fI -r„) 



iu n + fi - t(k) - S(iw„) 



(4) 
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where Y>(iuj n ) is the converged cluster self-energy, k is any 
vector in the original Brillouin zone and fi, v label cluster 
sites. All quantities plotted in this paper are derived 
from this lattice Green function. This is different from 
the periodization of the self-energy in Refs. 0, 0, 0] . 

Since spectra in and near a Mott insulator are expected 
to be highly singular, we used the exact diagonalization 
method [12| to solve the cluster-bath Hamiltonian (Eq.^J 
at zero temperature. This has the advantages of comput- 
ing dynamical quantities directly in real frequency and of 
treating the large U regime without difficulty. The self- 
consistency condition (Eq. |2J) was imposed on the imagi- 
nary frequency axis with a cutoff frequency larger than all 
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FIG. 1: Density of states JV(lj) at half-filling (n = 1) for 
several values of U. The dashed curves are N(u>) with AF 
long-range order. 



the bath energies. In the present study we used N c = 4 
sites for the cluster (minimum number of sites reflecting 
the full square lattice symmetry) and Nb = 8 sites for 
the bath. Four-fold symmetry of the cluster leads to four 
independent bath parameters. All energies are measured 
in units of t. Results are plotted with a broadening pa- 
rameter of O.li. 

In order to see how the Mott transition occurs in the 
presence of short-range correlations and how the differ- 
ences between short- and long-range correlation effects 
evolve with interaction strength, we present in Fig. ^ 
the density of states N(uj) as computed from the para- 
magnetic (P) (solid curves) and the antiferromagnetic 
(AF) (dashed curves) states at half-filling (n = I). The 
two states (P and AF) are selected by imposing appro- 
priate constraints on the bath parameters. When the 
correlation length is larger than the size of a cluster, 
long-range effects in the single-particle spectrum are ex- 
pected to be well described by the solution with the cor- 
responding long-range order. The difference between the 
two solutions at low energies is largest at weak coupling 
(U = 4t,5t) and decreases at strong coupling (U > 8t). 
In the weak coupling case, N(oj) has only a small dip at 
the Fermi energy compared to the AF solution, indicat- 
ing that short-range correlations available in a cluster of 
size N r = 4 are not long enough to lead to an insulating 
gap [Iil,[l4(- A s ^ i s increased to 6t, two sharp bands 
begin to develop at the gap edge and become well sep- 
arated from the Hubbard bands centered around ±Z7/2. 
At this strength of U, the two solutions (P and AF) start 
approaching each other at low energies. When U = 8t, 
an insulating gap is well defined between the low energy 
inner bands, and long-range correlations (dashed curve) 
sharpen these bands but they stay essentially at the same 
energy. Thus short-range correlations begin to dominate 
the low energy Physics at intermediate to strong cou- 
pling. At this value of U the low energy inner bands 
are clearly distinguished from the Hubbard bands com- 
ing from purely local correlations. Since these bands are 
missed in the single site DMFT ^| and smoothly evolve 
into sharp spin density wave-like bands when AF spin 
order is allowed, they are attributed to short-range spin 
correlations. When U is further increased beyond the 
bandwidth of 8t (not shown), the size of the insulating 
gap increases, but the relative shape and position of the 
inner and the Hubbard bands are unchanged. 

Next we investigate the evolution of a Mott-Hubbard 
insulator into a correlated metal upon a small electron 
doping at U = 8t. It is clear at half-filling (Fig. Efa)) 
that A(k,oj) (and also N(oj) in Fig. Q shows four sig- 
nificant features near (0,7r) and (tt/2,tt/2). The four 
peak structure in A(k,uj) was previously found by Moreo 
et al. I6J and Preuss et al. [17| in their high resolution 
Quantum Monte Carlo (QMC) results for the half-filled 
2D Hubbard model. As will be shown shortly, the ex- 
istence of the additional bands apart from the Hubbard 
ones turns out to be crucial for the systematic physical 
picture of the pseudogap phenomenon. 
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FIG. 2: Spectral function A(k,ui) along some symmetry directions in the paramagnetic (upper row) and the AF (lower row) 
states for several small dopings with U = 8t. The dashed lines represent the Fermi energy. 



The most striking feature of even a tiny doping into a 
Mott-Hubbard insulator (Fig. |2fb) and (c)) is that one 
of the inner bands at to > seems to disappear imme- 
diately while the Fermi energy jumps to the other band 
that persists at w < 0. Some of the weight is trans- 
ferred near the Fermi energy at w > 18]. Further- 
more, one already sees here the pseudogap phenomenon. 
Especially for n — 0.94, one clearly sees that the peak 
dispersing from low binding energy near (7r/2,7r/2) be- 
comes sharper close to the Fermi energy while near (n, 0) 
the corresponding peak becomes broader and less intense. 
The Fermi energy tends to sit towards a minimum in- 
stead of a maximum near that point. The suppression 
of weight near the Fermi energy eventually disappears 
at about 20% doping. The appearance of a low energy 
peak together with one at high binding energy is con- 
sistent with exact diagonalization results for the t — J 
model 0, ^| . We also show clear differences between 
the momentum resolved spectra of the paramagnetic and 
the AF states. The main differences for the AF case at 
half-filling (Fig.^d)) are that the low energy bands are 
more intense than in the paramagnetic state (Fig. G^a)) 
and that they bend back near (0, 7r) and near (7r/2, 7t/2), 
reflecting the underlying AF symmetry The results of 
Fig. |2fd) are strikingly similar to those obtained in a 10 
site V-CPT calculation |2Cj. The more drastic differences 
between the paramagnetic and AF solutions occur at fi- 
nite doping. While in both cases the Fermi level immedi- 
ately jumps to the closest low energy band upon doping, 
Figs, life) and (f) |2l| show that, in the AF solution, the 
four peak structure is preserved and the excitations at 
the Fermi energy are quite sharp. 
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FIG. 3: Density of states N(u>) for several doping levels with 
U = St. The solid, dotted, dashed, and long-dashed curves 
correspond to n = 0.96, 0.9, 0.84, and 0.74, respectively. The 
inset shows N(oj) near the Fermi energy. 



Fig ■ El shows the density of states N (oj) for several dop- 
ings. The pseudogap found here (inset) occurs over a 
scale 0.3i — OAt and vanishes for lar ge e nough doping 
in agreement with previous works |22l l23||. However, we 
have checked that at large U the scale of the pseudogap, 
as defined by the distance between the peaks closest to 
zero in Fig. |3J is almost independent of U in agreement 
with Ref. |24|. It does not scale as J, contrary to Ref. [2^| . 
Furthermore, the pseudogap size increases with decreas- 
ing doping, consistent with several recent experiments 
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FIG. 4: Imaginary part of the self-energy for (a) n — 0.84 
at k = (0, 52/64)/tt (solid line) and (30/64, 30/64)/tt (dashed 
line), and (b) n = 0.74 at k = (0, 49/64) /tt (solid line) and 
(29/64, 29/64) /tt (dashed line). X= (0,tt) and H= {tt/2, tt/2). 



in the cuprates Q, |2{|. The pseudogap feature disap- 
pears not because the two peaks coalesce but because of 
a gradual filling up of the spectral weight between the two 
peaks. Note the evolution with doping of the asymme- 
try about the Fermi energy. It is qualitatively consistent 
with what is expected for a doped Mott insulator [Isj . 
We also wish to point out clear evidence of correlation 
between the existence of pseudogap at the Fermi level 
(inset) and the separation of the low energy bands from 
the Hubbard ones. The appearance of a pseudogap near 
n = 0.8 found in the single-particle spectrum coincides 
with the downturn of the uniform cluster spin suscepti- 
bility (not shown) , suggesting that the pseudogap has to 
do with spin singlet formation of correlated electrons. 

The significant differences of spectral properties at 
n = 0.84 and n = 0.74 (Fig. can be understood in 
terms of self-energies [24| • Fig. Ufa) shows the imaginary 
part of the self-energy for n — 0.84 at the Fermi wave vec- 
tors j2|| near (0,7r) and (7r/2,7r/2). The scattering rate 
becomes large near X in particular, leading to the pseu- 
dogap in the single-particle spectral function. However, 
Fig. 01 b) shows much smaller self-energy in magnitude 
for n — 0.74, leading towards the Fermi liquid at large 
doping. Clearly, upon doping a Mott insulator into a 
metallic state, quasiparticles disappear since the absolute 
value of the imaginary part of the self-energy is a local 
maximum at the Fermi energy instead of minimum as in 
a Fermi liquid picture (see caveat The large scat- 

tering rate leading to the pseudogap found here is rem- 
iniscent of that due to the presence of the large density 
of states for scattering at the van Hove singularity in the 
weak coupling case. With a finite second-neighbor hop- 
ping t' , the crossing of the non- interacting Fermi surface 
and the AF Brillouin zone boundary leads to additional 
scattering processes 0. We now discuss this case. 

We study the asymmetry of hole- and electron-doped 
cuprates in their spectral function in the t — t' — U model 
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FIG. 5: Spectral function A(k, u) along some symmetry di- 
rections in the t — t' — U model with t' jt — —0.3 for 5% (a) 
hole- and (b) electron-doping. The dashed lines represent the 
Fermi energy. 



with t'/t = —0.3 [23 for 5% hole- and electron-doping. 
As in the case of t' = 0, the pseudogap shown in Fig. [S] 
appears in the low energy band caused by short-range 
spin correlations, which is well separated from the Hub- 
bard bands persisting at high energies. However, the 
locations where the pseudogap occurs are opposite for 
hole- and electron-doping. For hole-doping, A(k,u>) is 
gapped near k = (tt, 0), while for electron-doping it is 
gapped near k = (tt/2, tt/2), in agreement with ARPES 
experiments [2!J, y(J, IU Li3 • This is clearer on the cor- 
responding density plots in Fig. that show the asym- 
metry of hole- and electron-doping at u) — consistent 
with ARPES experiments. The combination of the above 
two plots indicate that the pseudogap is associated not 
only with the suppression of the spectral weight, but also 
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FIG. 6: Density plot A(k, 0) in the t - t' - U model with 
U = 8t, t'/t — —0.3 for 5% hole- (left) and electron-doping 
(right), i.e. n = 0.95 and n = 1.05, respectively. The current 
plots were obtained by taking the average of A(k, u) in a 
(—0.24, Q.2t) window around the Fermi energy. 

with its shift from the Fermi level. Since the spectral 
function A(k,u>) is obtained on a small cluster in the 
paramagnetic state, short-range spin correlations suffice 
for the pseudogap at both dopings. This is in agreement 
with CPT results I24| and exact diagonalizations on the 
t — t' — t" — J model |33j , but in disagreement with recent 
DCA results [34] where AF long-range order (or a large 
AF correlation length) is needed for a pseudogap under 
electron doping. 

It is remarkable that including only the nearest- 
neighbor correlations into this two-dimensional model is 
enough to capture some of the details of the spectral 
properties that are missed in the single site DMFT. We 
stress that at weak coupling and in two dimensions, a dif- 
ferent mechanism involving large AF correlation lengths 
causes the pseudogap 0, 0> HE EH] ■ That mechanism, 
relevant for electron-doped cuprates near optimal dop- 
ing, is clearly not involved in the present strong coupling 
calculations on small clusters. Although the sizes of clus- 
ters accessible to computation limit the resolution of the 
spectral properties, we expect the above picture for the 
strong coupling pseudogap to be robust with respect to 
an increase in the cluster size. This is confirmed by our 
recent low (but finite) temperature CDMFT+QMC cal- 
culations |37| where at intermediate to strong coupling a 
2x2 cluster accounts for more than 95% of the correla- 
tion effect of the infinite size cluster in the single-particle 
spectrum. 

Our calculations were done at zero temperature, while 
the pseudogap behavior is observed at finite temperature 
and broken-symmetry states appear in most of the zero 
temperature phase diagram. Nonetheless, our present 
results would be most relevant when the temperature is 



smaller than the characteristic energy scale for the pseu- 
dogap or in regions where broken-symmctrics are absent 
or are destroyed at zero temperature. In fact the lat- 
ter situation w as p recisely realized in a series of exper- 
iments [38l |3|| ^3 011 the cuprates, which have shown 
a metal-insulator crossover near optimal doping by sup- 
pressing the superconducting phase with a high magnetic 
field, revealing the underlying pseu dog ap state at zero 
temperature. Recent experiments [4l[ also show that 
when superconductivity is suppressed by Ni doping, the 
pseudogap survives at very low temperature in the corre- 
sponding paramagnetic state. Our results could be also 
relevant even in broken-symmetry states, as long as the 
sizes of the corresponding gaps are much smaller than the 
characteristic energy scale for the pseudogap. Even if we 
work at T — 0, the limited size of our clusters prevents 
long-range order unless we add mean-fields. Hence our 
calculation mimics the normal state. 

In conclusion, the long standing problem of the evolu- 
tion of a Mott-Hubbard insulator into a correlated metal 
upon doping has been examined in the two dimensional 
Hubbard model by using the Cellular Dynamical Mean 
Field Theory (CDMFT) which incorporates short-range 
spatial correlations. At half- filling these correlations cre- 
ate two additional bands besides the familiar Hubbard 
ones. Even a tiny doping into a Mott-Hubbard insula- 
tor [t! — 0) causes the Fermi energy to jump to one of 
these bands and the spectral weight to be suppressed im- 
mediately in a k dependent way, such that a pseudogap 
appears near (tt, 0) while a peak survives near (ir/2, 7r/2). 
When AF long-range order is present at finite doping, the 
picture is quite different, demonstrating that at strong 
coupling the pseudogap in the particle-hole symmetric 
model is a short-range effect. Including particle-hole 
asymmetry with a second neighbor hopping t' gives spec- 
tral weights that are similar to those observed in hole- 
and electron-doped cuprates. This provides a systematic 
physical picture for the emergence of the strong coupling 
pseudogap phenomenon in doped Mott-Hubbard insula- 
tors and the consequent non-Fermi liquid behavior that 
arises solely from short-range spin correlations [22ll23ll2^ | 
and large scattering rates, without any symmetry break- 
ing. 
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